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ABSTRACT 


Recent visual satellite photographs of certain intense convective 
storms have revealed concentric wave patterns. These waves appear as 
undulations in the top of the cirrus outflow of such storms, and are 
usually centered about the area of the most intense convection. They 
appear in both mid-latitude and tropical storms. This report analyzes 
two mid-latitude and three tropical cases of these waves. 

A model for the generation and growth of these waves is proposed. 
The proposed initial generating mechanism is similar to the effect 
noticed when a pebble is dropped into a calm pond. The penetration of 
the tropopause by overshooting convection is analogous to the pebble's 
penetration of the water's surface. The model for wave growth involves 
instability due to the wind shear resulting from the cirrus outflow. 
This model is based on an equation for the waves ' phase speed which is 
similar to the Helmholtz equation. It, however, does not assume an 
incompressible atmosphere, but rather assumes density is a logarithmic 
function of height. 

Finally, the model is evaluated on the two mid-latitude and three 
tropical cases. The data indicate that shearing instability may be a 
significant factor in the appearance of these waves. 
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INTRODUCTION 


1.1 Scope and Purpose 

Since the advent of meteorological satellites, interest in severe 
storm cloud-tops has increased. This is probably because the cloud- 
tops are usually all that can be seen of a storm from a satellite. This 
limited view, however, has exposed a number of new and striking 
phenomena. This report examines yet another one of these newly 
revealed phenomena . 

In a limited number of intense convective storms, the top surface 
of the cirrus canopy exhibits a concentric wave pattern. These waves 
appear to be centered about the area of the most intense convection. 

They have been observed in both the mid-latitudes and the tropics. 

This report proposes a physical and mathematical model, which 
under similar environmental conditions, would produce this type of wave 
pattern. The model is then evaluated by analyzing several case studies, 
two from the mid-latitudes and three from the tropics. 

1.2 Historical Background 

Since these waves are a relatively new discovery, research 
concerning them is scarce. However, a great deal of study has been 
carried out on other cloud-top characteristics of severe storms. In 
addition, much has been written on billow clouds and other gravity shear 
wave-type phenomena, not necessarily associated with intense convection. 
Pertinent developments in these two areas will be discussed in the 
following sections. 

Severe Storm Cloud-Top Studies . In one of the earlier studies on 
storm tops. Roach (1967) used photographs from U-2 aircraft and other 
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ISOTHERMS (“O 

OUTLINE OF CLOUD 

SUGGESTED STREAMLINES OF FLOW WITHIN CLOUD 

Figure 1.1. Model for parcel motions within a quasi-steady state 
overshooting dome. (From Roach, 1967). 

synoptic data to analyze severe thunderstorm cases over Oklahoma. Using 
these data he deduced a model for air flow in an overshooting dome 
(Fig. 1.1). He also observed waves in many of the cases he studied. 
These waves, however, were on the domes themselves and thus, had much 
shorter wavelengths than those observed by satellites. He occasionally 
observed a warm trench, as much as 1 km deep, around the perimeter of 
an overshooting dome. This may be a significant factor in the genera- 
tion of longer concentric waves (see section 3.1). Lee (1971), in 
analyzing this same data, also found the warm trench to be a significant 
feature. He, however, noticed it in both severe and non-severe storms. 

Purdom (1971) analyzed NOAA-1 and ATS-3 satellite photographs of 
severe storm tops. He found a correlation between tornado occurrence 
and a pause in the cirrus anvil growth. Fujita (1972) studied ATS-3 
photographs of a severe thunderstorm complex over Kansas. He found a 
correlation between tornado location and the location of an overshooting 
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Figure 1.2. Apollo 9 photograph of waves atop the cirrus outflow of a 
thunderstorm. Photograph was taken over Columbia, South 
America, March 8, 1969. (From Fujita, 1972). 


dome. The relative height of the dome, however, did not appear related 
to the tornado's formation. One interesting aspect of this paper is 
that it contained an Apollo 9 photograph of an anvil top over South 
America (Fig. 1.2). This photograph clearly reveals waves concentric 
about the storm's center. There was, however, no further discussion of 


these waves. 
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Fujita (1973) proposed a model by which a tornado would form when 
an overshooting dome collapsed. This was based on satellite and air- 
craft observations and laboratory experiments. Fujita (1974) further 
investigated overshooting dome characteristics. He likened the over- 
shooting dome to the bulge that would appear on a water surface when an 
upward shooting jet of water was injected below the water's surface. 

In the course of these investigations, he never observed concentric 
wave patterns. In fact, he made special note of the uniformly smooth 
appearance of the cirrus tops in all cases. This was in spite of the 
presence of numerous overshooting domes. He did, however, describe 
the effects of the rapid collapse of an overshooting dome. These 
effects appear to be a feasible means of wave generation. This is 
elaborated on in section 3.1. 

Arn (1975) studied anvil growth rates from ATS-3 and SMS-1 
photographs of seven severe storms over the United States. Again 
concentric waves were not observed. 

Arnold (1975) discussed concentric wave patterns observed in DMSP 
photographs of tropical storms. He stated that these patterns had been 
observed in DMSP photographs for some time. In fact, satellite meteo- 
rologists at Guam had been quite successful in using the focus of these 
waves for locating storm centers. He also postulated that they were 
gravity shear-type waves. 

Gravity Shear Wave Theories . As was mentioned earlier, a great 
deal has been written about gravity shear waves. However, only those 
publications utilized in creating the model in section 3.2 of this 
report will be mentioned. 
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Haurwitz (1941) described the derivation of the "Helmholtz" wave 
equation for an incompressible atmosphere. This equation is also 
derived by using various other methods in Haltiner and Martin (1957), 
Haltiner (1971), and Gossard and Hooke (1975). 

Gossard and Hooke (1975) also discussed stability criteria for a 
number of models for gravity shear waves. Much of their discussion is 
also contained in Drazin and Howard (1966) . Einaudi and Lalas (1974) 
explored wave stability for Kelvln-Helmholtz waves. 

The model described in section 3.2 of this report was first proposed 
by Gossard and Munk (1954) and further elaborated on in Gossard and 
Hooke (1975). Neither publication, however, contained a discussion of 
stability implications for this particular model. 
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OBSERVED WAVE CHARACTERISTICS 


Figures 2.1 through 2.5 are visual satellite photographs of cirrus 
outflow layers resulting from intense convection. These photographs 
display wave patterns with the following characteristics: 

1. The waves are visible as undulations in the top of the cirrus 
outflow. 

2. The waves appear to form concentric arcs or rings. 

3. The focal point of these arcs or rings appears to be the center 
of the most intense convection. This is characterized by 
either an overshooting dome, or in the case of Tropical Storm 
Ellen, by the apparent circulation center. 

4. In the Nebraska and Kansas cases the waves appear in the 
direction of greatest cirrus outflow (Figs. 2.1 -2.2). 

5. Areas of apparent amplification may occur away from the center 
of convection (Fig. 2.3). 

6. The range of observed wavelengths is 3.5 km (Fig. 2.5) to 14 
km (Fig. 2.2). 

7. The wavelength is not necessarily constant for a given cirrus 
outflow region. In the Nebraska case (Fig. 2.2) the wave- 
length increases further out from the dome center. In 
Tropical Storm Ellen (Fig. 2.3) the wavelength appears 
smaller further out from the storm center. 
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Figure 2.1. SMS-2 visual photograph taken over Kansas at 23^8Z, 
June 18, 1975- 







Figure 2.3» Sky lab visual photograph of Tropical Storm Ellen (Pacific 

1973)- 
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Figure 2.4. DMSP visual photograph of Typhoon Bess taken at 
October 10, 1974. 
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3. 


THEORIES OF WAVE GENERATION AND GROWTH 


3.1 Wave Generation 

The concentric orientation of the waves suggests generation at a 
point source. This point source would be located at the focus of the 
waves. As was mentioned in section 2.0, the focus of the waves gener- 
ally coincides with the center of the most intense convection. In the 
mid-latitude storms, this is marked by an overshooting dome penetrating 
deep into the stratosphere. In the tropics, this region is also char- 
acterized by deep convective penetration. Here, however this may be 
associated with either the eye wall of a tropical storm (Figs. 2.3- 
2.4), or the overshooting dome of an individual cell (Fig. 2.5). In 
all cases, however, it appears that convection penetrating through the 
tropopause is the key to generating the waves. 

Assuming this, the problem now is to determine what facet of this 
deep convection causes the waves. One guess may be that the waves are 
similar to the ripples caused when a pebble is dropped into a calm 
pond. The rapid penetration of the dome through the tropopause is 
analogous to the pebble's penetration of the pond's surface. In such 
a case, the wavelength would be expected to be similar in scale to the 
diameter of the penetrating body. The generating mechanism, however, 
would only be present during the initial penetration. Once the dome 
had reached a quasi-steady state (vertical and horizontal extent remain 
constant), this means of wave generation would cease. Therefore, some 
type of Ins'-ablllty would have to exist in order to sustain wave 
activity after the dome had reached a quasi-steady state. The Nebraska 
case as well as the tropical storm cases Indicate continued wave 
generation even though the penetrative convection appears to be 
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Figure 3-1* SMS-2 visual photograph taken over Nebraska at 0102Z, 

June 19, 1975- 

Figures 2-5 and 3*1 exhibit waves that appear to have been generated by 
this mechanism. 

Another possible means by which the waves could be generated would 
be by the periodic collapse and re-growth of the dome itself. This 
periodic dome growth could continually set up waves by generating wake 




14 



Figure 3.2. Model for parcel motions caused by a collapsing turret, 
(From Fujita, 1974). 

cirrus similar to the "turret eater" described by Fujita (1974) (see 
Fig. 3.2). This model can also be described by the pebble in the pond 
concept, except that repeated tropopause penetration continues to 
generate waves. Further wave instability is not necessary in this 
model. One drawback, however, is that the concentric nature of the 
waves requires the periodic collapse and growth of a single dome in 
roughly the same place. There is no evidence of such stationary 
oscillating domes presented by Fujita (1974) or Shenk (1974). Even 
more obscure is the means by which such a regularly periodic dome 
oscillation could be generated and sustained. 

The last proposed mechanism for wave generation is by quasi-steady 
state dome. A diagram of the particle motions within such a dome is 
given in Figure 3.3. Notice that this model is very similar to that 
proposed by Roach (1967) (see Fig. 1.1), except that it includes a 
warm trench along the dome perimeter. Roach (1967) did observe such a 
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Warm 

Figure 3.3. Alternate model for motions within a quasi-steady state 
dome. 

trench around a quasi-steady state dome. However, he chose not to 
include it in his model diagram. 

The creation of this trench is most simply explained by the "parcel 
method." A parcel of air located near the overshooting dome top 
(Fig. 3.3) will be colder than its environment. This negative buoyancy 
will cause a downward acceleration of the parcel. As it moves downward, 
it will also accelerate horizontally outward from the dome. A local 
high pressure created by the cold overshooting dome causes this 
acceleration Figure 3.4 shows this type of motion on a single 
turret observed by Fujita (1974). 

Once the parcel sinks adiabatically below the level of neutral 
buoyancy, it becomes warmer than its environment. Thus, this downward 
"overshooting" along with the horizontal parcel motion causes the 
creation of the warm trench. The mechanism responsible for the 
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Figure 3.4. Observed cloud particle motions on an overshooting turret, 
(From Fujita, 1974). 

formation of the warm trench would also be responsible for the 
generation of the waves further out from the dome. In fact, the 
warm trench may be viewed as merely a section of such a wave. 

By the preceding description, the phase speed of the trench 
(relative to the dome) is always zero. The phase speed of the 
waves further out, however, does not appear to be zero (see case 
studies). Again, further wave instability is not required in this 
model . 

Of the three generating mechanisms just described, the first 
appears to be the most reasonable. As will be shown in the case 
studies, shearing instability does appear to exist for waves in the 
cirrus outflow region of the storm. Thus, wave motion will persist, 
even after the dome becomes quasi-steady state. 


/ 
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Cirrus 

Outflow 


Figure 3.5. Diagram of typical concentric wave pattern. 

Model . Before describing the details of the model, let us examine 
the characteristics of the observed phenomenon. Figure 3.5 shows the 
features observed in visual satellite photographs. There are two impor- 
tant features to be considered in Figure 3.5. First, notice that the 
wavelength changes farther out from the dome. To simplify the model we 
will assume a mean wavelength over a given Ax, where Ax is measured along 
a specified radial emanating from the center of the dome. Secondly, if 
the cirrus outflow is diverging in two dimensions, then the horizontal 
speed of the outflow, measured along a given radial, will be inversely 
proportional to the distance from the dome center. To simplify this, 
we will assume a mean outflow rate over Ax. Since we are dealing with 
cirrus outflow resulting from intense convection, we will assume the 
cirrus cloud top to be at the tropopause. Figure 3.6 contains the 
vertical profiles for 9, T, p, and u for the model. Notice the two- 
dimensional flow in the xz plane is assumed. 
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Figure 3.6. Vertical profiles of T, 9, u^, and p for stability model. 

Layer 1 is the layer of cirrus outflow. It is located immediately below 
the tropopause. Layer 2 is the layer immediately above the cirrus out- 
flow and the tropopause. 

The following additional assumptions will be applied to (3.5) 
through (3.9). 

1. = 0, and v = 0. This is a result of the earlier assumption 
that all the motion takes place in the xz plane. This 
eliminates (3.6) from the basic set of equations. 

2. The tropopause is horizontal for undisturbed motion and 

z = 0 describes the height of the tropopause. This implies 
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that there are no significant motions resulting from a slope in 
the tropopause. 

3. The top layer is unbounded. This imposes a boundary condition 
at z = “ which is discussed later. 

4. The bottom layer is not infinitely deep but there are no effects 
from the lower boundary. Haurwitz (1941) described this as- 
sumption and showed that if the depth of the lower layer is 
greater than .4X, less than 1% error is introduced by this as- 
sumption. 

5. w =0. This means that the mean vertical velocity relative 

o 

to the tropopause is zero. 

3po ^Po 

6. -r — = 0, -r — = 0. Since p is the average pressure on a hor- 

3x dt o 

izontal surface, the first of these statements follows by 

definition. The second merely implies that the pressure at 

a given level does not change with respect to time. 

^Po ^*^o 

7. r — = 0, T — = 0. Again, the first statement is true by defi- 

dX ot 

nition. The second states that local changes in p with re- 

o 

spect to time are negligible. 

9u 8u 

8. ^=0, ^ - 0. Again, the first statement results from the 
chosen method of averaging. The second implies that the average 
motion at any level does not change with respect to time. 

~ ® within each layer. This has already been described by 
the selected vertical profile for u^ in Figure 3.6. 

10. The effects from the earth's rotation are negligible (f = 0) . 
Using these assumptions, (3.5) through (3.9) are rewritten as: 
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/ 8u ^ . 3P n 

/9w 3w\ ^ « 

“o “o 3?j+ 1T + '’* ■ ° 

M+U |E.+ „^+p |iL+o |» = o 
3t o 9x 9z 9x o 9z 

l£+u =yRt/^+ u ^ + 

9t o 9x 9z ^ I 9t o 9x 9z / 


Rewriting (3 . 12 ) : 


9p . 9p . 9u 9w 

9t 9x ^ 9z *^o 9x 9z 


(3.10) 

(3.11) 

( 3 . 12 ) 
(3.13) 


(3.14) 


Assiome hydrostatic atmosphere: 

3p 

^.-Ppg (3.15) 

Substituting (3.14) and (3.15) into (3.13): 


■" "o It ■ <“o8» = Ij 


C»U , 

^ + p 


dW 

o 9z 


Define operator ^ as: 


5- = i- + 1_ 

Dt 9t '‘o 9x 


Since w = 0 and v = 0: 
o o 


(3.16) 


(3.17) 


( ^ , L- 2_ 

V*^^/o o 9x Dt 


Rewriting (3.11) using (3.17): 


Dw 1_ ^ p_ 

Dt Po 9z S = 0 (3.18) 


- i 


I 
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Taking ^ of equation (3.18): 

d!w + 1_5_ /i2.\+S_^ = 0 

Dt2 Po Dt \^z) Pq Dt 
Rewriting (3.13) using (3.15) and (3.17): 


(3.19) 


P 




_ Dp g_/ Dp \ ^ ^ „ 

^ Dt p^yRT^Dt/ ^ ^ 


YRT 


p 8z 


(3.20) 


Define: 


® \^YRT P^ 3z / 


(3.21) 


Gossard and Hooke show that N is actually the Brunt-Vaisala frequency 
squared. Also, by use of the first law of thermodynamics for adiabatic 
motion and the equation of state, it can be shown that: 


. 4 ^ 

0 dz 


Substitute (3.21) into (3.20): 


^ Dp = 22. + N^w 

Po Dt p yRT Dt 


Substitute (3.23) into (3.19): 


(3.22) 


(3.23) 


or: 




(3.24) 


Employing the techniques of Gossard and Hooke (1975), define the follow- 
ing parameters: 
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U E 


W E 


P E 



(3.25) 

(3.26) 

(3.27) 


where : 


= density at some reference level 
-6z 

Po = e Ps 
6 = some constant 

These substitutions are a direct result of the density profile illustrated 
in Figure 3.9. That is, varies exponentially with respect to height. 
(See section 6.0 for further discussion on the errors associated with 
this assumption.) This transformation is made to maintain constant co- 
efficients in our equations while still allowing compressibility. 

Using (3.25), (3.26), and (3.27), we will now develop a new set of 
basic equations based on (3.10), (3.16), and (3.24). First, substitute 
(3.24) and (3.27) into (3.10): 



or: 


DU + ^ 

Dt pg 3x 


0 


(3.28) 


/O 

- v-'* 


OC.\ J 

aiiu 


/O 0-7N 
J 


, - / O O \ . 

lUU } • 




Dt 


F/ 1 


+ ) P 

\2po 

3z 

YRT 9z / 


= 0 


(3.29) 
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Again,- following Gossard and Hooke (1975), define: 


r = 


1 

1 ° + 


= i (_a_ _ 

2 \yRT g / 


2po 3z yRT 
Substitute (3.30) into (3.29): 

Substitute (3.25), (3.26), and (3.27) into (3.16) 

I 2F + ™+ ('3 B 

p yRT Dt 3x \3z yRT 2p 3z / 


Substitute (3.30) into (3.32): 


1 DP ^ 3U . 

+ -T T 


P^YRT Dt 3x 




(3.30) 


(3.31) 


(3.32) 


(3.33) 


For convenience, the three basic equations just discussed are rewritten 
below: 


0 

Dt p 3x 
s 


W = 0 




(3.28) 

(3.33) 

(3.31) 


Utilizing these three equations, we will now investigate the possible 
forms of solutions for U, W, and P for each layer of the model. A 
subscript of "1" will denote the lower layer and a subscript of "2" 
the upper layer. 

Let us first investigate solutions for the upper layer. Assume 
solutions of the form: 
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U2 = [y(x-ct)] (3.3A) 

^2 ~ ^ 2 ^ [y(x-^t)] (3.35) 

?2 = 026 ^2^ cos [yCx-ct)] (3.36) 

To limit the possible solutions, two boundary conditions were Imposed. 
First, at X = 0 and t = 0, let W2 = 0. Second, at z = “, let W2 = 0. 

A solution set for the lower layer may be obtained in a similar fashion 
resulting in: 

cos Cv(x-ct)] (3.37) 

sin [p(x-ct)] (3.38) 

cos [y(x-ct)] (3.39) 


The second boundary condition does not necessarily apply in the lower 
layer since it is not infinitely deep. Let us look at w^^ as a function 




(3i+|6i)z 

w^ = C^e sin [y(x-ct)] 


Notice that an amplification in w^ as we move downward from z = 0 
(tropopause) will occur only if 3^^ < ~ ^ ^1’ fora of amplification 

in the vertical of the lower layer should be minimized as much as pos- 
sible. Otherwise, if our wave is effecting a significant Az in the 
lower layer, large amplitude errors in the vertical will be introduced. 
Nevertheless, the condition 6^^ ^ ” Y *^1» require 3^^ > 0 for 

the lower layer. That is, for certain 3j^ < 0, amplification of w^^ will 
still not occur as we move downward from z = 0. Thus, 3j^ can be positive 
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or negative so long as the magnitude of a given negative 3 ^ is not 
extremely large , 

Now, directing our attention to the upper layer only, substitute 
(3.34), (3.35), and (3.36) into (3.28), (3.33), and (3.31): 


I 

S <'=-“02> - paj °2 ■ “ 



(3.40) 

1 

yRI^ ”2-*2-<=2( 


0 

(3.41) 



■ (r 2 - 62 ) 

(c-Uo2)I>2 = 0 

(3.42) 

Solve 

(3.40) for D 2 and substitute into 

(3.41): 



YRT 2 j ^2 ^2 ^2 ( 


0 

(3.43) 


or: 


X2^2~^2 


(») 


= 0 


(3.44) 


where : 


^2 yRT, 


(C-U02) -1 


Solve (3.44) for A 2 and substitute into (3.40) 

<=2 / »2+''2 \ , 


— D = 0 


° 2 ' Ps, “2 


Solve (3.46) for D 2 and substitute into (3.42) 

2 2 

2 9 9 ^2 ~^2 2 

U (c-u. ) - r- — (c-Uo,) = 0 


(3.45) 


(3.46) 


(3.47) 


Solve (3.47) for 82 ' 
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- 1 / ^2 1 
Xl - ^ 


(3.57) 


Solve (3.56) for and substitute into (3.52); 


( - M (c-uo,) C + — D = 

V Xj^y / ”11 psj^ 1 


= 0 


Solve (3.58) for and substitute into (3.54): 

2 , .2 2 ^^1 "^1 ^ , >2 - n 

y (c-Uqj^) - (c-Uqj^) - 0 


(3.58) 


(3.59) 


Solve (3.59) for 6^^; 


3i = ± 



r 2 1 


+ Xl 

N 2 

I 

2 " ^ 

- (c-^Uoj^) 

I 


) 


(3.60) 


Using (3.56) and (3.58) to substitute for and in terms of in 
(3.37) and (3.39), the solution set for the lower layer may now be 
written as: - 


/r -8 \ 8 z 

Ui = — j e cos [y(x-ct)] 


3,2 


Wi = e sin [y(x-ct)] 


( r \ 3 z 

e ^ cos [y(x-ct)] 


(3.61) 

(3.62) 

(3.63) 


Utilizing the boundary condition for a surface of discontinuity described 
by Haurwitz (1941), the following equations for each layer are obtained: 


3(Pj^-P2) 3(P3^-P2) 3(Poi+Po2) 

— H — + “°i — 5 ;^ — “1 — 55 — ’'“1 


3(P;^-P2) 


3y 


+ V, 


3 (poj^-po2) 


3y 


+ W, 


3(P^-P2) 

°1 dz 


+ W, 


3(PO]^-Po2) 


(3.64) 


3z 


= 0 
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aCP^-P^) 9(P^-P2) 3 (Poi-Po2) 9(Pi-P2> 

+ u r + U„ r + V 


3t 


O 2 9x 2 8x 


°2 


+ V, 


9(Po]^"^Po2^ 9 (Pj^ P 2 ) 9 (Poj^~Po2^ 


+ w 


+ W„ 


O 2 9z "2 3z 


= 0 


2 3y 

Employing the assumptions discussed earlier, (3.64) and (3.65) 


3(P^-P2> 3<Poi-Po2> 

+ U r + W 


3t 


o^ 3x ’ "1 3z 


= 0 


9(P^ 9(Po,-Pq2) 

3t ^°2 9x ^^2 3z “ ° 

From (3.51) and (3.63): 

P^ - P^ = c cos [y(x-ct)] 

where: 

/r -6 \ (t z /s,+r,\-s, 

0 = c;p,, <c-u„,) z - C'P,, (c-u„^> ■ 

From the hydrostatic equation, (3.15): 

Pn ~P^ ' 8^ ) 

°1 °2 °1 °2 


We will now choose p and p such that: 

Si S 2 

p = p measured at the upper surface of the tropopause 

"1 °1 

p = p measured at the lower surface of the tropopause 
®2 °2 

Substitute (3.71) and (3.72) into (3.70): 

Po - ?o = - 82 (Ps "^s ^ 

12 1 ®2 

Substitute (3.68) and (3.73) into (3.66) and (3.67): 


3iZ 


I 

(c-u )o-c'e g(p - p^ ) = 0 
°1 1 "1 ®2 


(3.65) 

reduce to: 

(3.66) 

(3.67) 

(3.68) 
z 

(3.69) 

(3.70) 

(3.71) 

(3.72) 

(3.73) 


(3.74) 



31 


-3„z 

,(c-u )o-C'e ^ g(p -p ) = 0 
°2 ®1 2 


(3.75) 


Substitute (3.69) back into (3.74) and (3.75) and rearrange terms: 


-3-z 


(c-u ) (c-u ) p 
°1 


°2 ®2 






(3.76) 


3,z . .2 

C'e ^ ^ ^s 

°2 ®2 






(c-u_ ) (c 




~^s ^ 
®2 


Equate (3.76) and (3.77): 


(c-u ) (c-u )p 
°1 °2 ®2 


(¥) 


(3.77) 


(c-u )^p ( ^ "“s ^ 

°2 ® 2 \ ^2 / ®1 ®2 




or: 


(3.78) 


(c-u )p- (c-u )p^ +g(P -P^) =0 (3.79) 

°1 \ / °2 ®2 \ ^2 / ®2 


Using the equation of state and Poisson's equation, we will now 

substitute for p and p in terms of 0. and 0 (each measured 
®1 ®2 12 

immediately below and above the tropopause respectively). First, from 
the equation of state: 

(3.80) 


s, RT 
1 s. 


(Subscript s^^ indicates lower surface of tropopause.) 
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(Subscript S 2 indicates upper surface of tropopause.) 
Secondly, from Poisson's equation: 



where: 

K = constant 

Substitute (3.82), (3.83), and (3.84) into (3.80) and (3.81): 



Again, since p 


P 

®1 ®s. 


p 



(3.81) 


(3.82) 


(3.83) 


(3.84) 


(3.85) 


(3.86) 


(3.8/) 


(3.88) 
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where : 


K’ = 


= constant 


RK RK 

Substitute (3.87) and (3-88) into (3.79): 


(3.89) 


(c-u )‘ 
°1 


0 + (c-u ) 6 

S2V Xi / 02 


T 2+ e 2 


+ g (0 -0 ) 

®2 


= 0 (3.90) 


We will now define a mean horizontal wind speed between the two 
layers : 


u +u 
°1 °2 

% = 2 


(3.91) 


Align the coordinate system such that: 


u = 0 
o 


Therefore: 


u = -u 
°2 °1 


(3.92) 


and : 


2u = Au 

Of 


(3.93) 


where Au is the difference between the mean wind speeds of each 
o 

layer, measured across the tropopause. Thus, under this new coordinate 
system, (3.90) becomes: 


2 /^1 ^ l \ 2 

(c-u ) 0 + (c+u ) 0 

°1 ^2 \ H / °1 ®1 


r +8 

— - + g (0 - 0 ) = 0 (3.94) 

X2 S2 


where : 


^2 i ^2 ^2 


N. 


- y 


(c+u )‘ 
°1 


'2 


(3.95) 
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(3.96) 


(3.97) 


(3.98) 


and 3^ and are defined by (3.60) and (3.57). 

Equation (3.94) will be the basic equation used to determine wave 

stability. As was mentioned earlier, the real component of is 

always positive. However, no such restriction was placed on 6^. 

Therefore, for a given u there are two possible solutions for c, 

°1 

(assuming a solution does exist for both $^^>0 and 3j^<0) . One 
feature, however, must be maintained, and that is, c must be a con- 
tinuous function of u . Appendix B discusses the selection of 

°1 

3j^>0 or 3j^<0. 

Wave Stability . The ultimate goal of this model is to accurately 


depict wave Instability with respect to time. From equations (3.38) 
and (3.35), it can be seen that if c is purely real, no amplification 
will occur, and thus, the waves wi‘11 be neutrally stable. If, however, 
c is complex, growth will occur based on the following: 

Let : 

c = c + c . i 
r 1 

Sxnce : 

W^, “ sin [y(x-ct)] (3.99) 

and: 

c , c . , u , X and t are all real 
r 1 
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Therefore: 

[(|ix-yc^t) - (pc^t)i] (3.100) 

or: 

W^, oc sin (m + ni) (3.101) 

where : 

m = yx - yc^t (3.102) 

n = -yc^t (3.103) 

From Kreyszlg (1968) : 

sin (m + ni) = sin m cosh n + i cos m sinh n (3.104) 


If and in (3.38) and (3.35) are assumed to be purely real, the 

imaginary term on the right side of (3.104) is not physically signif- 
icant. Therefore, the amplification rate of W 2 and is cosh 

(-yc^t), where c^ is the imaginary component of c. For convenience, 
we will define the instability amplification factor, ij;: 

fp = cosh (-yc^t) (3.105) 

Figure 3.7 contains traces of c , |c,|, and 8, versus Au . 

The characteristics displayed in this example are also found in each 
of the case studies. Those characteristics considered most significant 
are: 

1. In the c versus Au curve, c is a continuous function 

r o r 

of Au . 
o 

2. In the |c^| versus Au^ curve, |c.| is greater than zero 
only over a range of 4u^. Also, |c^| displays an additional 
point of neutral stability (|c.| = 0) within this range. 

This point occurs at c = . 



Cr(m /s 
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3. In the 3^^ versus Au^ curve, 3^^ does go from a positive to 

^Uq 

a negative value. At c = ~^» 3j^ = ± “>• This is difficult 

to interpret physically as far as 3^^ itself is concerned. 

However, it does not introduce discontinuities in c. and 

1 

c^ since the coefficient of 3^^ is zero at this point. 

Thus, we see that the model can be used to infer instability regimes. 


This will be elaborated on further in the case studies. 



4. DATA ANALYSIS TECHNIQUES- 
4.1 Mid-Latitude 

The data used to analyze the mid-latitude cases (Nebraska and 
Kansas) is summarized in Table 4.1. Figure 4.1 shows the area of 
coverage for the upper • air, radar, and surface observations. 


Table 4.1. Mid-Latitude Data Sources 


Type of Data 


Time of Observations Data Source 


Satellite Photographs 
(IR and Visible) 


Once every 15 to 30 min 
from 12Z June 18, 1975 
to 12Z June 19, 1975 


SMS-1 and SMS-2 photo- 
graphs received at the 
Satellite Field Service 
Station (SFSS) in Kansas 
City, Missouri 


Upper Air Soundings 


Radar Reports and 
Surface Observations 


NMC Facsimile Maps 


12Z June 18, 1975 and 
OOZ June 19, 1975 


Once every hour plus 
specials from 12Z 
June 18, 1975 to 12Z 
June 19, 1975 

In accordance with — 
NAMFAX schedule from 
12Z June 17, 1975 to 
12Z June 20, 1975 


Computer plotted sound- 
ings received at the 
Severe Local Storms 
Forecast Center in 
Kansas City, Missouri 

Teletype reports re- 
ceived at F.^ E. Warren 
Air Force Base, Wyoming 


Facsimile products re- 
ceived via NAMFAX cir- 
cuit GF-10201 at the 
Department of Atmo- 
spheric Science, 
Colorado State 
University 


Scale factors for the satellite photographs were derived from map 
overlays such as the one shown in Figure 4.2. These overlays were added 
to the prints made by the SFSS in Kansas City. They were not on the 
original negatives. The following example illustrates two distance 
measurements on a blow up of this photograph. 
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Figure 4.1. Map showing area of data coverage for the mid-latitude 
cases . 
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19JN75 32A-H 





# 
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I h 


Figure 4.2, SMS-2 visual photograph with map grid, taken at 0102Z, 
June 19, 1975. 
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In Figure 4.2, calculate the east-west map scale factor, E , 

s 

for the Nebraska area. To do this, measure the length of the 

southern border of the Nebraska panhandle on the map. Ax , 

Pm 

and compare this to its actual length on earth. Ax : 

Pa 

Ax =27.5 mm 
Pm 

Ax = 165 km 


Ax 


m 


E = A - C 

s Ax t 

p 6x10 


Again in Figure 4.2, calculate the north-south map scale 

factor, N , for Nebraska by similar- means, this time using 
s 

the western border of the panhandle as a reference: 

Ay = 26 mm 
Pm 


Ay = 222 km 

Pa 




N = 


■ m 


® ^^p 8.5x10^ 


Measure the distance between two numbers on the legend of the 


original print (AN^ in Fig. 4.2) and between the same two 
numbers on the blow-up (AN in Fig. 4.3). Compare these 

D 

distances to calculate an enlargement factor, F: 

48 mm 


F = 


22 mm 


= 2.18 


Multiply E and N by F to get the corresponding scale 
s s 

factors for the blow-up, E_ and N_ : 

D D 



2 



Figure 4.3. Blow up of Figure 4.2 without map grid. 
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Eg = (2.18) X 


6x10^ 2.75x10^ 


Ng = (2.18) X 


8.5x10^ 3.9x10^ 


Measure the east-west diameter of the cloud on the blow-up. 

Ax , and multiply it by to obtain the cloud's actual 

east-west diameter. Ax : 

c 

a 

Ax = 135 mm 
*^B 

Ax = (135 mm) x (2.75x10^) 
c 

a 

= 371 km 

Measure the north-south diameter of the cloud on the blow-up. 

Ay , and multiply it by to obtain the cloud's actual 
^B ^B 

north-south diameter. Ay : 

a 

Ay = 83 mm 
^B 


Ay = (83 mm) x (3.9 x 10 ) 
c 

a 


= 324 km 


The data used to calculate wave stability criteria, according to 
the model is section 3.2, were taken from the nearest upper air sounding 
(time and space) to each storm. The following is an example of how 
these data are incorporated into the computer program which solves 
equation (3.94) for c in terms of Au^ and X: 

1. Read the temperature and pressure immediately above and below 
the tropopause (T^, P 2 and T^, in Fig. 4.4). 

2. Read the temperature and pressure of the next reported level 

above and below the tropopause (T , P„ and T , P 

np Up low low 


in Fig. 4.4). 
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6CX) 

700 

800 

900 

1000 

1050 


Figure 4.4. Temperature sounding from North Platte, Nebraska taken 
at OOOOZ, June 19, 1975. 


Calculate a mean temperature in each layer using: 


Ti+Ti 


"l = 


low 


T +T 
2 2 


T = 
2 


_up_ 


Using the hydrostatic equation, calculate Az^, between 


and P^ , and between P^ and P^ 

low up 


Using the 

equation of state to substitute for p in the hydrostatic 


equation, the solutions for Az^ and become: 


RT^ 

A-1 =-7 


In 


low 
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5. Using Poisson's equation, calculate 9 for each of the four 
pressure levels. For example: 


6 . 



Solve for and N 2 using the finite difference form 

of equation (3.22): 




Once these parameters are determined, equation (3.94) is solved 

for c, for selected Au and A. This is done by iteration using 

o 

"Muller's Method with Deflation," which is contained in Colorado State 

University IMSL subroutine ZANLYT. 

As was described in section 3.2, the real component of c is the 

phase speed of the wave relative to the coordinate system where u^ = 

0. Therefore, the phase speed of the wave relative to the dome is: 

c* = c + u * 
r o 


where: 


u* + u* 


* = 


and c*, u*, u* , and u* are now measured relative to the dome. 


O 0-, 


Therefore: 


.* = 


u* + u* 
°1 °2 


+ c 


2 
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Also, Au^ is the same regardless of the coordinate system. Thus: 


Au = u 
o o. 


- u 


= u* - 


Au 

. ‘ . c* = — 4- u* + c (4.1) 

2 O 2 r 

This equation is used to calculate c*, the phase speed relative to 

the dome, when Au , u* , and c are known. Since u* is the 

o O 2 r . °2 

stratospheric wind speed relative to the dome, it is calculated by 
subtracting the dome's motion from the observed stratospheric wind on 
the nearest sounding. 


4.2 Tropical 

The data used to analyze the tropical cases is summarized in 
Table 4.2. The method of analysis is essentially the same as for the 
mid-latitude cases. However, the map scale factor for DMSP photographs 
is given as 1:1.5x10^ and is assumed uniform in all directions through- 
out each photograph. This is permitted because the satellite is more 
directly overhead in the DMSP photographs than in the mid-latitude 


SMS-2 photographs. 
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Table 4.2. Tropical Data Sources 


Type of Data 

Time of Observations 

Date Source 

Satellite Photographs 

Once per case, 
usually near 2300Z 

DMSP photographs 
received at 
observing sites in 
the western Pacific 

Upper air soundings 

OOZ and 12Z daily 

"Daily Weather Maps 
with Synoptic Data 
Tabulation," pub- 
lished by Japan 
Meteorological 
Agency , Tokyo , 

Japan 

Daily Synoptic Maps 

OOZ and 12Z for surface 
analysis, 12Z for upper 
air analysis 

"Daily Weather Maps 
with Synoptic Data 
Tabulation," pub- 
lished by Japan 
Meteorological 
Agency , Tokyo , 

Japan 



5. 


CASE STUDIES 


5.1 Nebraska 

Between 2300Z on June 18, 1975 and 0230Z on June 19, 1975, a severe 
thunderstorm cell moved through Central Nebraska, spawning eight tor- 
nadoes over a nearly continuous path 170 km long (Fig. 5.1). This storm 
is the subject of our first case study and will henceforth be referred 
to as the "Nebraska" storm. 

Figure 5.2 contains SMS-2 visual photographs showing the early 
evolution of the storm. The rapid growth of this storm is evident in 
these photographs. The radar reports from Grand Island, Nebraska during 
the most active period of the storm are plotted in Figure 5.3. Notice 



Figure 5.1. Path of tornadoes spawned by Nebraska storm on 
June 18, 1975. 
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2000Z 


2030Z 



2100Z 


2200Z 


2130Z 



2230Z 



Figure 5.2. SMS-2 visual photographs taken over the Kansas-Nebraska 
area on June 18, 1975. 
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• 
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• 
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(Only the tornado cell was 
reported on this special) 



0230Z 19 Jun 75 


Figure 5.3. Grand Island, Nebraska radar summaries during period of 
tornado activity. 


that the tornado producing cell was in the southeast section of the 
echo region and that its top remained between 57,000 ft and 61,000 ft. 
Also, note that the cell movement was southerly at approximately 15 kts 


throughout the period. 
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The concentric waves atop the cirrus outflow region of this storm 
first appeared on the 2317Z SMS-2 photograph (Fig. 5.4). They remained 
visible until 0147Z (Figs. 5.5 -5.7), when sunset prevented the taking 
of any further visual photographs. Table 5.1 summarizes the wavelengths. 



Table 

5.1. Nebraska 

Case Cloud 

Measurements 


QUADRANT 

TIME 

WAVELENGTH 

AMPLITUDE^ 

RADIUS^ 

OUTFLOW RATE^ 



(km) 

(km) 

(km) 

(m/sec) 

North 

2317Z 

— 

— 

70 

18 


0102Z 

— 

— 

185 

8.5 


0132Z 

— 

— 

200 



0147Z 

— 

— 

— 


East 

2317Z 

10 

1-3 

75 

23 


0102Z 

11 

1.5-3 

220 

22 


0132Z 

11-14.5 

.35-. 85 

260 

28 


0147Z 

11-14.5 

.25-. 40 

285 


South 

2317Z 

— 

— 

30 

11 


0102Z 

— 

— 

100 



0132Z 

— 

— 

— 



0147Z 

— 

— 

— 


West 

2317Z 

10 

— 

60 

17.5 


0102Z 

— 

— 

170 

5.5 


0132Z 

— 

— 

180 

11 


0147Z 

— 

— 

190 



See Appendix A. 

2 

Relative to southeast dome. 

3 

(A Radius) t (A time). This is the radial speed of the cloud 
boundary. 






Figure 5.4. SMS-2 visual photograph taken over Nebraska at 2317Z, 
June 18, 1975. 


I 
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Figure 5.5. SMS-2 visual photograph taken over Nebraska at 0102Z 
June 19, 1975. 



Figure 5.6. SMS-2 visual photograph taken over Nebraska at 0132Z, 
June 19, 1975. 


Figure 5.7. SMS-2 visual photograph taken over Nebraska at 0147Z, 
June 19, 1975. 
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amplitudes, cirrus radii, and computed cirrus outflow rates for each 
photograph. 

The vertical profiles of wind and temperature from the OOOOZ North 
Platte, Nebraska sounding on this date, are given in Figures 5.8 and 
5.9. Table 5.2 lists the parameters from this sounding that were used 
to solve equation (3.9A). 



Table 

5.2. 

Thermal 

Parameters for 

Nebraska Case 

^2 

= 212°K 


= 125 mb 



up 


up 




^2 

= 212°K 

"2 

= 150 mb 

6^ = 364. 7°K 

2 -4 -1 

N 2 = 4.51x10 sec 

"l 

= 206. 7°K 


= 153 mb 

6^ = 353. 6“K 

N^^ = 2.12xl0"''^sec"^ 


= 215. 3°K 


= 200 mb 




low low 


Before discussing the results of equation (3.94), let us first 
examine some of the important features exhibited by the data just 
presented : 

1. The radar reports (Fig. 5.3) indicate that the sounding was 
taken less than 50 km north of the storm's center. The 
satellite photographs also indicate that the sounding very 
likely passed through the cirrus outflow layer of the northern 
portion of the storm. 

2. The observed wind maximum near the tropopause (Fig. 5.9) 
apppears to be the result of the intense cirrus outflow at 
that level. The shallowness of the wind maximum and the fact 
that it is near the layer where cirrus outflow would be 
expected, help support this conclusion. The presence of a 





Height { km ) 



Temperature (°C) 


Figure 5.8. Temperature sounding from North Platte, Nebraska taken 
at OOOOZ, June 19, 1975. 



u ( m/sec) V (m /sec ) 

Figure 5.9. Wind profiles from North Platte, Nebraska sounding at 
OOOOZ, June 19, 1975. 
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jet stream to the west and north of the area (Fig. 5.10), 

however, indicates that this could also be the result of a 

larger scale feature. In any event, Au across the tropopause 

o 

is nearly 30 m/sec in this area. 

3. The eastern quadrant of the cirrus outflow contained the 

most distinct wave pattern throughout the period (Figs. 5.4- 
5.7). It also had the strongest outflow rate (Table 5.1). 




) 


Figure 5.10. Winds at 35,000 ft over the Western United States at 
OOOOZ, June 19, 1975. 
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4. Since the radius measurements were made with respect to the 
southeastern dome (Table 5.1), the southern and eastern outflow 
rates are most likely good indications of the radial wind 
speed at cloud level. In the western and northern quadrants, 
however, outflow from new cells (Figs. 5. 3-5. 7) enhances the 
motion of the cloud boundary. Wind speeds at cloud level 
cannot, therefore, be accurately inferred by cloud expansion 
rates in these quadrants. 

5. The stratospheric wind relative to the dome (Figs. 5.3 and 

5.9) is approximately westerly at 10 m/sec. Therefore, the 

average Au at the eastern cloud edge is 15 m/sec. 
o 

6. The top of the southeastern dome, according to the radar 
reports (Fig. 5.3), remained between 2 km and 5 km above the 
tropopause during the period for which waves were visible. 

The dome's diameter, estimated from the 0132Z satellite 
photograph (Fig. 5.6), was approximately 11 km. 

7. The wavelengths in the eastern quadrant remained fairly 
constant throughout the period (Table 5.1), averaging about 
11 km. In the 0132Z and 0147Z photographs (Figs. 5.6 -5.7), 
however, waves in the far eastern portion of the outflow 
increased to a wavelength of about 14.5 km. 

8. The amplitude of the waves showed a marked decrease from 
0102Z to 0147Z (Table 5.1). 

The thrust of all this is that there appears to be sufficient 
evidence that (1) a substantial tropopause penetration did occur in 
the form of an overshooting dome, (2) the cirrus outflow created some 
degree of wind shear across the tropopause, (3) the diameter of the 
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dome is of the same scale as the wavelength, and (4) the waves occurred 

in the area of apparent maximum wind shear across the tropopause. 

Using the data from Table 5.2, we will now evaluate equation (3.94) 

for various X and Au . In doing so, we will attempt to evaluate 

o 

the model's applicability to the case. Figure 5.11 shows the variation 

of |c^| (magnitude of the imaginary component of the phase speed) and 

c* (real component of the phase speed measured relative to the dome) 

with respect to Au , when X is held constant. Notice that each 

o 

wavelength is unstable only over a certain range of Au for which a 

o 

purely real solution for c does not exist. For values of Au either 

o 

side of this range, neutral stability prevails. This result is quite 
unlike that obtained by the incompressible Helmholtz equation 
(Fig. 5.12), where |c^| continues to increase once a critical Au^ 
is attained. 

From section 3.2 we saw that the stability amplification factor, 

\p, equals cosh (yc^t). Figure 5.13 contains a sample plot of ip vs 

X for a given Au^ and t = 30 minutes. Notice that the plot has 

two humps. This is consistent with the two hump configuration of |c^| 

vs Au in Figure 5.11. Plots of the right hump only and the left 
o 

hump only of the \p vs X curves for a series of Au^ and t = 30 

minutes, are given in Figures 5.14 and 5.15. 

Using Figures 5.14 and 5.15, curves of Au vs X for a constant 

o 

ip are plotted in Figure 5.16. This graph shows that for a given Au , 

o 

a certain range of X's will be the most unstable. The two humps 
again appear, this time as two axes of maximum instability. 

Let us now return to the observations to evaluate the model's 


results. First of all, the average wavelength is 11 km. Secondly, 
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Figure 5.15. Left hump of ip versus A 
(Nebraska case). 
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Figure 5.16. Au versus A curves for selected \p at t = 30 min (Nebraska case). 
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based on observed cirrus outflow rates and the sounding data, Au in 

o 

the area of the waves appears to be between 15 m/sec and 30 m/sec 

(Table 5.1 and Fig. 5.9). Figure 5.16 indicates that 11 km is one of 

two most unstable wavelengths when Au = 22 m/sec and when Au = 

o o 

56 m/sec. Thus, for the lower predicted value of Au , at least, the 

o 

data does appear to be in fair agreement with the model. However, if, 

in fact, Au = 22 m/sec. Figure 5.16 would indicate that X = 3.7 km 
o 

would be even more unstable than A = 11 km. The apparent absence of 

this wavelength may be related to the fact that the dome's diameter is 

approximately 11 km. A dome of this size would likely produce higher 

amplitude waves at A = 11 km than at A = 3.7 km. Thus, even though 

ip is slightly larger for A = 3.7 km, the 11 km wavelength waves still 

prevail since their initial amplitude was much larger. 

There are two remaining observed wave characteristics that the 

model does not appear to predict. The first of these is the apparent 

drop in wave amplitude from 0102Z to 0147Z. The second is the increase 

in A as the waves progress further eastward (Figs. 5.6 and 5.7). 

These two features may be related. Let us assume that near 0102Z the 

shear across the tropopause decreased rapidly to Au^ = 10 m/sec. 

Figure 5.16 would indicate that the most unstable wavelength would now 

be approximately 3.5 km and wavelengths above 6.5 km would be stable. 

If, however, this decrease in Au occurred rapidly, iL for A = 3.5 km 

o 

would still be relatively small since {p increases with respect to 
time. Thus, we will attribute the absence of visible A = 3.5 km waves 
to the lack of time for amplification. The increase in A further east 
may also be explained by this rapid drop of Au^ to 10 m/sec. 

Figure 5.17 shows the dispersion curve for Au^ = 10 m/sec. Notice 
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that the longer wavelengths move at a greater phase speed than the 
shorter wavelengths. Thus, just as in the case of a sea swell, the 
longer wavelengths will move out ahead of the shorter wavelengths. 

This appears to be the case in Figures 5.6 and 5.7. 

The conclusion, therefore, for this particular case, is that the 
stability model of section 3.2 can be used to describe some of the 
observed wave characteristics. However, lack of precise measurements 
over a substantial length of time prevent us from making a more 
rigorous evaluation of the model. This is especially true with respect 
to the features described in the preceding paragraph. 

5.2 Kansas 

On the same date as the Nebraska storm, an intense squall line 
also developed over western Kansas and Oklahoma. We will refer to this 
storm as the "Kansas" case. This system produced several funnel clouds 
in its northern section, but none touched down. In the southern 
section, near Guyman, Oklahoma, golf ball sized hail was reported at 
about 0140Z on June 19, 1975. 

Again, the early evolution of this storm system can be seen in 
Figure 5.2. The radar reports from Garden City, Kansas are given in 
Figure 5.18. Like the Nebraska storm, this system developed very 
rapidly. The average cell movement throughout the period was south- 
southwesterly at approximately 45 kts. The cell apparently responsible 
for the waves in this case was located near the Kansas-Oklahoma border 
on the 0030Z radar summary. Presumably, this cell was also responsible 
for the hail near Guyman, Oklahoma. 

Wave patterns on the cirrus tops of this system appeared most 
evidently on the 2348Z SMS-2 photograph (Fig. 5.20). Photographs taken 
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after this revealed very little wave activity (Figs. 5.20 - 5.24). The 
next photograph before 2348Z was taken at 2317Z (Fig. 5.19). It too 
contains no distinct wave features. Table 5.3 summarizes the key 
measurements taken from Figures 5.19-5.24. 

The upper air sounding from Dodge City, Kansas, taken at OOOOZ on 
this date is given in Figures 5.25 - 5.26. Table 5.4 lists the 
parameters from this sounding that were used to solve equation (3.94). 

The important features to be noted in the data for this case are: 

1. The radar reports (Fig. 5.18) indicate that the sounding was 
taken about 110 km east of the squall line's axis. The 
satellite data, however, indicates that the sounding passed 
just east of the cirrus outflow layer. 

2. There was no dramatic wind shear observed at the tropopause 
level of the Dodge City sounding (Fig. 5.26). 

3. The area of wave activity was also the area of maximum cirrus 
outflow (Fig. 5.19). 

4. Only outflow rates perpendicular to the squall line were 
computed. 

5. The stratospheric wind relative to the cells (according to 
the radar reports and the Dodge City sounding) was approxi- 
mately northerly at 8 m/sec (Figs. 5.18 and 5.26). Therefore, 

the average Au at the eastern cloud edge was 15 m/sec (this 
o 

measurement was made at the maximum eastward cirrus outflow). 

6. The diameter of the dome at the focus of the waves was 
approximately 8 km (Fig. 5.19). 

As in the Nebraska case, there again appears to be sufficient 
evidence that (1) a substantial tropopause penetration occurred. 
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Figure 5.19. SMS-2 visual photograph taken over Kansas at 2317Z, 
June 18, 1975. 
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Figure 5.20. SMS-2 visual photograph taken over Kansas at 2348Z, 
June 18, 1975. 




Figure 5.21, SMS-2 visual photograph taken over Kansas at 0002Z, 
June 19, 1975. 




Figure 5.22. SMS-2 visual photograph taken over Kansas at 0017Z, 
June 19, 1975. 
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Figure 5.23. SMS-2 visual photograph taken over Kansas at 0032Z, 
June 19, 1975. 






Figure 5.24. SMS-2 visual phoi 
June 19, 1975. 


taken over Kansas at 0047Z 
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Table 5.3. Kansas Case Cloud Measurements 


QUADRANT 

TIME WAVELENGTH 

AMPLITUDE 

RADIUS^ 

OUTFLOW RATE^ 


(km) 

(km) 

(km) 

(m/ sec) 

East 

2317Z 

— 

95 

14 


2348Z 5 

.5-1.2 

120 






17 


0002Z 

— 

135 

17 


001 7Z 

— 

150 

17 


0032Z 

— 

165 

17 


0047Z 

— 

180 


West 

2317Z 

— 

65 

8.5 


2348Z 

— 

80 

5.5 


0002Z 

— 

85 

11 


0017Z 

— 

95 

5.5 


0032Z 

— 

100 

5.5 


0047Z 

— 

105 


^Measured relative to the dome 

which appears 

to be at the 

focus of 

the waves 

in Figure 5.19. 




2 

(A Radius) 

T (A Time) . 





(2) cirrus outflow created some degree of wind shear across the 
tropopause, (3) the dome's diameter is of the same scale as the wave- 
length, and (4) the waves occurred in the area of apparent maximum 
Au^ across the tropopause. 

Using the data in Table 5.4, equation (3.94) was solved. Figures 

5.27 - 5.29 were constructed from these results. Again, the two humps 

in the \p versus X curves for constant Au result in the two axes 

o 

of maximum instability in Figure 5.29. Utilizing these two axes 



Height { km ) 



Figure 5.25. Temperature sounding from Dodge City, Kansas taken at 
OOOOZ, June 19, 1975. 



u(m/sec) v(m/sec) 


Figure 5.26. Wind profiles from Dodge City, Kansas sounding at OOOOZ 
June 19, 1975. 
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Table 5.4. Thermal Parameters for Kansas Case 


^2 

= 209. 0°K 

?2 

= 125 

mb 




up 


up 






^2 

= 208. 0°K 

"2 

= 150 

mb 

02 = 357. 9°K 

" 2 ' 

= 1.44x10 ^sec ^ 

"l 

= 207. 1°K 

"l 

= 154 

mb 

6^ = 353. 6°K 

^ 1 ' 

= 5.01xl0~‘^sec"^ 

"l. 

= 217. 9°K 


= 200 

mb 





low low 


one finds that 5 km is one of two most unstable wavelengths at 

Au = 13 m/sec and Au = 25 m/sec. This lower predicted value is 
o o 

again very close to the value of Au (15 m/sec) calculated from the 

o 

cloud motions and the sounding data. 

When Au = 13 m/sec on Figure 5.29, the lower axis of maximum 
o 

instability predicts that X - 2.3 km will also be a relatively un- 
stable wavelength. Again, the absence of visible waves of this wave- 
length may be due to their initially low amplitude relative to that of 
the 5 km wavelength waves. 

Finally, the waves were no longer visible after 2348Z. This was 

in spite of the apparent increase of Au^ to 17 m/ sec throughout 

the next hour. Figure 5.29 would predict most unstable wavelengths 

of 6.5 km and 3.2 km for Au = 17 m/sec. There are two possible 

o 

rationalizations for the absence of visible waves during this hour. 

The first is that the surge of outflow that caused the waves was 
not evident at the cloud edge until later. This means that a lag may 
exist from the time a surge first occurs in the inner portions of the 
cloud to the time it appears at the cloud edge. In other words, 
the Au^ in the area of the waves may have decreased already but this 
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Figure 5.27. Right hump of ip versus A curves for selected Au 's (Kansas case). 


= 30 min . 
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Figure 5.28. Left hump of ip versus X curves for selected Au 's (Kansas case). 





Figure 5.29. Au^ versus A curves for selected ip at t = 30 min (Kansas case) 


84 


will not be noticed at the cirrus edge until some time later. The 
second rationalization is that the Initial amplitude of the other 
predicted wavelength waves was so small that the amplification has not 
yet been sufficient to make them visible in the satellite photographs. 

5.3 Tropical Cases 

As was mentioned earlier, concentric waves have occasionally been 
observed in DMSP photographs of tropical storms for several years. 

Three cases from the western Pacific are shown in Figures 5.30- 5.32. 
These cases were selected from 156 visual DMSP photographs of Western 
Pacific tropical storms and typhoons. Of these 156 photographs, taken 
from June 1971 to January 1972, sixteen displayed some type of con- 
centric wave pattern. The three cases discussed in this section were 
chosen because of their proximity to an upper air sounding at the time 
of the photograph. The data pertaining to these storms are summarized 
in Table 5.5. Since only one photograph exists for each case, cirrus 
outflow rates could not be computed. Table 5.6 summarizes the thermal 
parameters from the closest sounding to each storm. 

The important features to be noted in the data for these cases are; 

1. Only tropical storm "kit" had a sounding within 600 km of 

Table 5.5. Observations from Tropical Cases 


STORM NAME 

TIME AND DATE 
OF PHOTOGRAPH 

WAVELENGTH 

(km) 

DISTANCE FROM 
CLOSEST SOUNDING 

Gilda 

2300Z June 24, 1971 

5-6 

950 km 

Rose 

2200Z August 12, 1971 

4.5 

650 km 

Kit 

2230Z January 4, 1972 

3.5 

less than 100 km 
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Figure 5.32. DMSP visual photograph of Typhoon Kit taken at 2230Z, 
January 4, 1972. 
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Table 5.6. Thermal Parameters for Tropical Cases 


Gilda 


T = 191. 9°K 

^2 = 

100 mb 




up 

up 





T^ = 192. 7°K 

"2 = 

112 mb 

02 = 370. 7°K 

CM 

CM 

= 4.99x10 ^sec ^ 

T^ = 203. 1°K 


150 mb 

0^ = 360. 4°K 


= 1.79xl0~^sec"^ 

T = 218. 5°K 

low 

^1, 

low 

200 mb 




Rose 






T = 194. 7°K 

^2 = 

100 mb 




up 

up 





= 193. 9°K 

^2 = 

104 mb 

02 = 376. 2°K 


-4 -1 

= 4.9x10 sec 

T^ = 203. 9°K 

"l = 

150 mb 

0^ = 370. 4°K 


= 2.5x10 ^sec ^ 

T = 219. 5°K 

low 


200 mb 




Kit 






T_ = 185. 9°K* 

^2 = 

80 mb * 




up 

up 





T^ = 185. 9°K 

"2 = 

100 mb 

02 = 359. 2°K 


= 5.15x10 ^sec ^ 

T^ = 204. 9°K 

"l = 

150 mb 

0^ = 355. 5°K 


= 3.5x10 ^sec ^ 

^ T = 219. 7°K 

low 


200 mb 





*Estimated values. The sounding ended at 100 mb yet did not report a 
tropopause. 
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its center. This sounding, however, did not report a 
tropopause. It too then, is not a very reliable source for 
the thermal parameters used in equation (3.94). 

2. No measurements of Au^ across the tropopause in each storm 
were available. 

In spite of these shortcomings, there is still evidence that 
(1) the focus of the waves in each photograph is characterized by 
penetrative convection and (2) the intense cirrus outflow would result 
in some wind shear near the tropopause in each case. 

Figures 5.33 - 5.35 contain graphs of the stability regimes 
calculated for each of these cases. In spite of the poor data base for 
each of the curves, they do display features consistent with the mid- 
latitude cases. First of all, the lower values for Au for which 

o 

each observed wavelength is one of two most unstable wavelengths, are 

quite close to those observed in the mid-latitude cases. Specifically, 

these values range from Au^ = 10 m/sec for Tropical Storm Rose to 

Au = 15 m/ sec for Tropical Storm Gilda. The curve for Tropical Storm 
o 

Kit predicts A = 3.5 km (observed wavelength) to be one of two most 

unstable wavelengths at Au = 11 m/sec. 

o 

The magnitude of ijj, predicted in each of these cases, however, 
varies considerably from case to case. For Tropical Storm Kit ip = 

1.035 when the observed wavelength is one of two predicted most unstable 
wavelengths. For Tropical Storm Rose, this value jumps to ip = 16 
where the observed wavelength is one of two most unstable. This large 
difference in ip from case to case is probably a result of errant 
data. As will be shown in section 6.0, the magnitude of ip is very 


sensitive to data errors. 



91 



in o o 


(UJM)\ 


Figure 5.34. Au versus A curves for selected ijj at t = 30 min (Typhoon Rose). 
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Unfortunately, the predicted values of Au^ in each case cannot 
be verified by observations. However, the mere fact that they are 
so close to those observed in the mid-latitude cases lends credence 


to the stability model. 



6. ERROR ANALYSIS 


6.1 Measurement Errors 

The data presented in this paper is, for the most part, subject 
to a great deal of uncertainty. This is especially true of the sounding 
data in the tropical cases. Even in the mid-latitude cases, the 
soundings were taken outside the area of visible wave activity. For 
this reason, a proposed range of expected values would be more or 
less a guess. Also, the distance measurements made on the satellite 
photographs are subject to error. These distance errors become even 
more significant when calculating amplitudes based on shadow lengths 
(see Appendix A). Rather than performing an arbitrary error analysis 
of each case study, we will use the Nebraska case data to demonstrate 
the sensitivity of the stability model to a selected amount of error. 

The Nebraska case was chosen because of its apparently superior data. 

We will also discuss amplitude errors resulting from inaccurate 
measurements and the perhaps invalid assumptions inherent in 
equation (A. 9). 

In the case of errors for the stability model, we will choose an 

error range of ± 1°C for each temperature measurement, T^ » T , T , 

■^iow ^ ^ 

and T^ • It can be shown that the stability characteristics would be 
up 

effected most if these errors were combined in the following two ways: 
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Error sets #1 and #2 applied to the Nebraska case at A = 11 km and 

t = 30 min result in curves #1 and #2 respectively in Figure 6.1. 

Obviously, the amplitude of ip for this case is greatly affected by 

these errors. The values of Au for maximum ip, however, are not 

o 

changed so drastically. For the lower hump in Figure 6.1, maximum \p 

is found at Au = 14 m/sec in curve #2 and at Au = 25 m/sec in 
o o 

curve #1. These values are still within the range of possible values 

in section 5.1. The upper hump shows that the maximum \p is found 

at nearly the same value of Au for both cases. The conclusion from 

o 

this is that the ranges of Au for which a given wavelength will be 

o 

unstable can be reasonably predicted using conventional data in the 
model. However, the degree of instability is far too sensitive to 
observational errors to permit accurate calculations of amplitude 
growth rates. 

The most serious amplitude errors are probably a result of the 
assumptions which are implied in equation (A. 9). Some of the suspected 
error producing assumptions are: 

1. The cirrus top is smooth except for the wave pattern. 

(Figure 1.2 shows that this is not necessarily the case.) 

2. The wave pattern is a sine wave. 

3. There is no refraction of diffraction of the sun's rays. 

4. The cloud's surface is totally opaque to the sun's visible 
radiation. 

The seriousness of these assumptions, however, cannot be accurately 
evaluated without a case which compares calculated values to direct 
amplitude measurements. The effects of measurement errors, regardless 
of these assumptions, are shown by the range of amplitudes given in 
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Tables 5.1 and 5.3. These ranges were produced by assuming measurement 
errors of ±1 km for both wavelength and shadow length. Obviously, only 
an order of magnitude estimate of amplitude is possible with these 
data. 

6.2 Stability Model Errors 

There were a multitude of assumptions made in the derivation of 
equation (3.94). Rather than attempting to evaluate each of these, we 
will look at the one assumption that separates this model from the 
incompressible model. That assumption is that is an exponential 

function of height rather than a constant. This is expressed in the 
following equation: 



Let us now see how this reduces the error in the lower layer of 

the Nebraska case. First, for the incompressible case, p for the 

o 

entire layer is estimated by p^ at the tropopause, (153 mb) . Using 

2 

the equation of state, this p^ is calculated to be .258 kg/m . If 

2 

we now go down approximately 2 km to 200 mb, p^ becomes .324 kg/m . 

Thus the error introduced at this level is 20 percent. On the other 

hand, the exponential model results in virtually no error in p^ at 

-4 -1 

this level, assuming 6 = 1.23x10 m in the lower layer. Slight 

errors will be introduced as lapse rates change in each layer, but 
these errors are negligible compared to those introduced by the 
incompressibility assumption. 

Before concluding our discussion of the errors in the stability 

model, we should mention one apparent inconsistency. This inconsistency 

occurs in the area of 3, = ± °°, which results at c = u 

1 o^ 


This is 
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similar to the infinity resulting at the critical level of other 
gravity wave models. Its physical significance is still in question 


at this point. 


7. 


SUMMARY AND CONCLUSIONS 


This report has examined concentric wave phenomena atop the cirrus 
level of intense convective storms. These waves appear in satellite 
photographs of both mid-latitude and tropical storms. They are usually 
centered about an area of deep penetrative convection. 

A model for the generation and growth of these waves was derived. 
This model proposes that the rapid penetration of the tropopause by 
overshooting convection initially generates the waves. Once generated, 
the waves are then further amplified by shearing instability. The 
instability equation was derived by perturbation techniques similar to 
those used to derive the Helmholtz equation. Incompressibility, how- 
ever, was not assumed. 

The approach used in creating the model for these waves had two 
basic steps. The first step was to determine which features occur in 
all of the cases observed. These Included penetrative convection, 
intense cirrus outflow, and a cirrus cloud top location near the tropo- 
pause. Utilizing these characteristics, a model containing these same 
features was developed. Once developed, the results of the model were 
then tested for consistency with respect to the observed conditions for 
each case. 

Two mid-latitude and three tropical cases were analyzed. The 
mid-latitude cases were observed from SMS-2, and thus several short 
time interval photographs were available. This, along with supporting 
radar and upper air data permitted calculations of the relative cirrus 
outflow rates which could then be used to approximate values for the 
wind shear across the tropopause. These calculated values were then 
compared to those predicted by the stability model. In this manner. 
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the validity of the model in the mid-latitude cases was evaluated 
quantitatively. 

The tropical cases, on the other hand, contained only one 
photograph of each case. Thus, values for the wind shear across the 
tropopause could not be calculated. Therefore, only a qualitative 
evaluation of the stability model was made in the tropical cases. 

The results from the case studies led to the following conclusions: 

1. These concentric waves are Initially the result of a rapid 
penetration of the tropopause by a deep convective storm. 

2. The cirrus outflow from the storm creates a wind shear across 
the tropopause. 

3. This wind shear is responsible for further wave growth due to 
shearing instability. 

Any conclusions regarding the validity of the wave stability model 
are relatively premature at this point. However, there are several 
distinct characteristics of the model that do warrant consideration: 

1. Incompressibility is not assumed and thus, errors created by 
this assumption in the Helmholtz equation are greatly reduced. 

2. The results of the model predict only a certain range of 
unstable wavelengths for a given wind shear across the tropo- 
pause. This is unlike the incompressible Helmholtz model 
which predicts increasing instability for all wavelengths 
shorter than the "critical wavelength." It would appear 
from observations, that the former results better describe 
atmospheric gravity shear wave phenomena. 
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3. The range of solutions obtained from the stability model in 

the case studies agrees favorably with observations, even when 
the data may be subject to substantial errors. 

This paper presents only a "first guess" as to the cause of these 
waves. In light of current observational data, it is an entirely 
feasible explanation. Only further study with more direct measurements 
of environmental conditions and wave motions can verify or disprove the 


model conclusively. 



8. SUGGESTIONS FOR FURTHER STUDY 


Obviously, a great deal of further testing of the stability model 
is needed before its true validity can be determined. In particular, 
better sounding data within the area of wave activity is needed. Also 
time lapse photographs would greatly aid in the verification of 
predicted phase speeds. 

Another important study of these waves would be a statistical 
analysis of their relationship to severe weather occurrences. This 
would require a much larger and more random sample than that presented 
in this paper. However, the few cases analyzed so far do show 
encouraging signs of some type of correlation. 

Finally, if the stability model is, in fact, valid, there is a 
possibility that these wave patterns could be used to deduce storm 
motions. This would be done by determining the storm moting necessary 
to create the observed wave pattern. Of course, temperature and wind 
profiles of the storm area would be necessary to make such a calcula- 
tion. Still, it is possible that even composite-type values of these 
parameters would be sufficient to deduce at least a direction of 
movement, and maybe even a range of possible speeds. 
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APPENDIX A 

AMPLITUDE MEASUREMENT TECHNIQUE 

Figure A.l contains a cross-sectional view of an idealized wave 
pattern. Figure A. 2 contains an overhead view of this same pattern. 

This idealized picture is used to derive a method for infering amplitude 
based on the shadow length in each wave. 

In Figure A.l and A. 2, the following parameters can be measured or 
obtained from solar tables, if the time and date of the photograph are 
known: 

length of the shaded portion of the wave 
x^ length of the bright portion of the wave 

X wavelength 

a 90° - (sun zenith angle) 

b 90° - (sun elevation angle) 

In Figure A. 2, we note that the wavelength as seen from the sun is 
A'. This can be calculated by using: 

X’ = — 

cos b 

Now, assuming that the sun's rays are straight parallel lines, we 
may write the equation for the ray passing through A and B as: 

z = mx + z^ (A.l) 

Where: 

m = slope of line AB 

Also the wave equation, as seen from the sun is: 
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z = z„ cos TTX (A. 2) 

^ A 

At A, the slope of the wave equation is the same as the slope of line 
AB. From (A. 2), this slope is: 

^2 ^ ~ ^ (A. 3) 

We now substitute (A. 3) into (A.l) at point A: 


Zj^ = - (tan a) + z^ 


(A. 4) 


Solve for z^ in equation (A. 2) at point A: 


2ir 


'l ~ ^2 


(A. 5) 


Substitute (A. 5) into (A. 4): 


2tt 

Z^ = cos ^ 


(A.6) 


Similarly at point B: 


z. = -(x, + X ) tan a + z„ 
4 1 s 3 


(A. 7) 


and : 


z, = 


cos 


[f^ '*1 ’‘s'] 


(A. 8) 


Therefore: 


2ir 


^ (xt + X )1 + X tan a 
LA -L S j a 


(A. 9) 


where from (A. 3) : 




(A. 10) 
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Equations (A. 9) and (A. 10) can now be solved for which is the 
amplitude of the wave (see Fig. A.l). 
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APPENDIX B 

SELECTION OF SIGN FOR 


The selection of 3^^<0 was, at first, determined by the 

process of illimination. Equation (3.94) was solved in each case, 
using both positive and negative values for In doing so, a certain 

trend became apparent. That is, for a selected sign for and certain 
ranges for Au^, equation (3.94) failed to converge on a solution. Thus, 
the sign of for which equation (3.94) did converge, within this same 
range of Au^> would be used. This sign would then be maintained until 
either g^ or its coefficient became zero. The points where g^ changes 
sign are illustrated in Figure 3.7. 

As a means of illustrating why equation (3.94) fails to converge 
for certain g^ and Au^ a typical case will be described. The solution 
of equation (3.94) for this case will be determined at u^^ = 0. 

At this point equation (3.94) simplifies to: 


~ ^(0S2^2 ^S2^l) + (0S]^r2 + 0s2ri)J + g (0s2~®si) = 0 (B.l) 


The following typical parameter values will be used: 


= 350° K 

®1 

Tj^ = 200° K 

2 -3 -2 

N^ = 2x10 sec 


0S2 = 370° K 

T 2 = 210° K 

N? = 3.8xl0"^sec”^ 


2 2 

The values of N^ and N 2 are both indicative of very thermally stable 
on scale analysis of this case, (B.l) becomes: 


2 

^ (0si32 + 0S2^i) 8 ^®S2 “ = 0 


(B.2) 


where : 



(B.3) 



2 2 

Notice in this case N 2 > . This is generally true since the air 

immediately above the tropopause is usually more stable than that immed- 
iately below. 

At this point, one further simplifying assumption is made. Again, 
based on scale analysis: 



(B.5) 



(B.6) 


The validity of these assumptions is not readily apparent at this point. 
However, the data from the case studies does indicate that solutions for 
c at Uqj^ = 0 are changed less than 10% under assumptions (B.5) and (B.6). 
Also, the results regarding the sign of 3^^, without these assumptions, 
are the same as those we will obtain utilizing them. Thus, (B.3) and 
(B.4) become: 



Now a wavelength of the same magnitude 


(B.7) 

(B.8) 

of those typically observed is 


chosen: 
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Therefore : 

2 -7 -2 

P = 3.95x10 m 

Under these circumstances, possible solutions for 3^ > 0 will be 

investigated. First, we will look for a real solution for c. Since 
2 2 

^2 ^ ^1 ’ minimum real solution is: 



Using the given values for N 2 and p , this minimum real c is: 
c = 71 m/sec 


Upon substituting this value into (B.7), (B.8), and then (B.2) it is 


found that the first term of the left. 




(0Sj^B 2 + 0S2^l) 


, is greater 


than the second term, g (0s2"®s-j^) • Also, any further increase in c 
causes even larger differences between these two terms. Therefore, a 
real solution for c when uqj^ = 0 and > 0 does not exists. 

Next the possibility of either a complex or purely imaginary solu- 
tion for c when 3^ > 0 will be investigated. Rewrite (B.2) under the 


conditions (B.5) and (B.6): 


0S'- 


^ 2 ^1 
- -T- 


(2 ^2 

r - j 


2 . 


(®S2 ®S]^) 


or: 


Tn 

L '“2' 




2 2 2 1 

O I ~ /Q _Q N 

*^'2 •' J ® '•''S2 


/t» N 
\a.-3] 


Since g(0s2~®Sj^) is real, then c must be the complex conjugate of the 
factor between the brackets. For example, let: 
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c = c + c.i 
r i 


(B.IO) 


where c and c, are real. Therefore: 
r i 


,22 1 22 2 -^ 

0S2(y c ) + 0si(y c^-N2 ) = + c^i 


(B.ll) 


Also : 


(Cr+CiD (c^-c^i) = g(0s2-esi) 

or: 

= g(0s2-0si) (B.12) 

2 2 2 

Now substitute for g, 0gj^, 9s2» » ^2 ’ from the given case 

into (B.ll) and (B.12): 

[(7.6xl0"^)c^-(7.3xl0“^)] + [(7.4xl0^)c^-(3.7xl0^)] ^ = 

(B.13) 

Cr + c^i 

2 

.'. = ± (196 - c^^) (B.14) 

Substitute (B.14) into (B.IO): 

2 ^ 

c = + (196 - ) + c^i (B.15) 

Substitute (B.15) into (B.13): 

{[ 7 . 6x10"^ ( 196-c^^) -7.3x10"^] + (7 .6xl0“^)c^i 1 + 

/ o % n _/■ I 

] [7.4x10 ° (196-c^ ) -3.74x10 ] + (74x10 ®)cii \ = 

2 

(196-c^^) -c^i 


(B.16) 
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Since and are real, from (B.14) the following condition is 
imposed : 

|c^| <14 (B.17) 

First we will investigate the possible solutions where c^ = 0 
(i.e., c is imaginary). 

From (B . 14) : 

c^ = ± 14 

If c^ = 14 and c^ = -14 are substituted into (B.16), a complex solution 
is obtained on the left side of the equation and a purely imaginary 
solution on the right side. Therefore, an imaginary solution for c in 
this case does not exist. 

The only remaining possible solution is where c is complex. To 
examine this possibility, only the sign of each term will be dealt with. 
That is , if a term is greater than zero , it will be represented by 
{POS}, and if a term is less than zero by {NEG}. 

First, we will examine the terms between the braces on the left 
side of (B.16). Notice, under condition (B.17): 

(7.4x10"^) (196-c^^) < 3.74x10“^ (B.18) 

-ft 9 ^ 

(7.6x10 ) (196-c^ ) < 7.3x10 ^ (B.19) 

Thus, the signs of (B=16) are: 

For c^ > 0: 
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[{NEG} + {POS}i]^ + [{NEG} + {POS}i]^ = 

(B.20) 

{POS or NEG} - {POS}i 
For < 0: 

[{NEG} + {NEG}i]"^ + [{NEG} + {NEG}l]^ = 

r (B.2i; 

{POS or NEG} - {NEG}l 


The formula for complex square root is defined as: 


J = m + ni 


v'J 




J I +m + (sign i) 




-m 


(B.22) 

(B.23) 


Notice, only with those roots that have a positive real component are 
dealt with. Therefore, (B.20) and (B.24) reduce to: 


[{NEG} + {P0S}i]^ + [{NEG} + {P0S}i]^ = {POS} + {NEG}i 


(B.24) 


[{NEG} + {NEG}i]^ [{NEG} + {NEG}i]'^ = {POS} + {P0S}i 


(B.25) 


Apply the form of solution given by (B.22) and (B.23) to equations 
(B.24 and (B.25) : 

{POS} + {P0S}i ^ {POS} + {NEG}i (B.26) 

{POS} + {NEG}i ^ {POS} + {P0S}i (B.27) 

Therefore, there are no complex solutions for c when Uq^ = 0 and 3^^ > 0 
for this case. 

The conclusion, from all of this, is that at Uq^ = 0, 3j^ > 0 will 
not necessarily have a solution. Therefore, 3^ 0 is selected to 
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determine the solution for c. Note that for 3^^ < 0, (B.26) and (B.27) 
do not necessarily result due to the difference in sign between 3^ 
and • 

As was mentioned earlier, c must be a continuous function of u^^. 
Since we have chosen 3^^ (real) < o at u^^ = 0, we must continue to use 
3^ (real) < 0 until either 3^^ (real) = 0 or the terms containing 3j^ 
as a factor have a real component equal to zero. The first condition. 



(B.28) 

(B.29) 


The latter condition is met at; 


c = Ug^ (B.30) 

As is demonstrated in Figure 3.5, 3^^ does shift sign when (B.28) is 
true, and again, where (B.30) is true. 
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